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Abstract Gauss' law and the equation of continuity must be satisfied in all 
materials be they solids, liquids or gases. Most materials are classified as simple 
materials; i.e., their electrical properties are linear, isotropic and homogeneous. 
Charge transport in these simple materials should be described by a constitutive 
equation known as Ohm's law. When Ohm's law is combined with Gauss' law 
and the equation of continuity, a differential equation for volume charge density 
relaxation results. The usual solution to this equation shows that charge decays 
exponentially with a relaxation time given by the material's permittivity divided 
by its electrical conductivity. Experiments show that good conductors follow 
this exponential decay but that poor conductors (insulators) tend to follow a de- 
cay that initially is more hyperbolic than exponential. This suggests that either 
Ohm's law is not valid for insulator materials or that a deeper understanding 
of Ohm's law is needed to explain charge decay in these less than good con- 
ductors. This paper examines the latter approach and shows that, when all the 
free charges within a simple material are taken into account, a new closed- form 
equation for unipolar charge decay can be derived which is valid for any simple 
material be it classified as a very good conductor, or a very poor conductor or 
anywhere in between. For good conductors the equation reduces to the standard 
exponential law of decay. For very poor conductors the equation reduces to the 
Vellenga-Klinkenberg modified hyperbolic law with the initial decay producing 
the characteristic Bustin hyperbolic law of decay. Explicit definitions for a good 
conductor and a good insulator are obtained and are used to define the range 
where explicit deviations from both of these hyperbolic laws occur. 

Introduction 

Electrostatics is the study of charges including their effects and how they af- 
fect the local environment. In other words, electrostatics examines the behavior 
of charges in various material media (solids, liquids and gases). In the ab- 
sence of magnetic influences, the relevant electrical properties of a material are 
its electrical conductivity a and its permittivity e given by e = eoCr, where 
eo = 8.85 X 10""'^^ F/m is the permittivity of free space and is the relative 
permittivity or dielectric constant of the material. 

When excessive charge builds up in a material, it can often have dire con- 
sequences such as initiation of a spark that results in a fire, explosion or other 
damage to the local environment. For example, in every operation from the sim- 
ple (e.g., walking) to the complicated (e.g., semiconductor manufacture) there 
is an operation time (e.g., walking step, ion-sputtering time, etc.) during which 
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a material (e.g., sole of shoe, semiconductor chip) can become charged. If the 
rate at which charge is deposited is greater than the rate at which charge can 
decay, a charge build-up on or within the material can occur. As a result, it 
is important to understand charge decay in sufficient detail so as to be able to 
predict if a material constitutes an electrostatic risk in a particular operation. 

In order to gain a more fundamental insight into charge decay, this paper first 
reviews the basic concept of applying Gauss' law and the equation of continuity 
to simple materials; i.e., materials that conform to Ohm's law. Then a more in- 
depth review of the physics which leads to Ohm's law is presented including the 
concept of electrical conductivity. This is followed by a more in-depth review 
of the partial differential equation that results when Ohm's law is combined 
with Gauss' law and the equation of continuity. After defining the concept 
of a material (solid, liquid or gas; conductor or insulator) as being electrically 
neutral in the intrinsic state, a unipolar charge species is introduced into a simple 
material. It is shown that this charge temporarily changes the overall electrical 
conductivity, a point either missed or ignored in past charge decay derivations. 
Once it is recognized that the conductivity has changed due to the charge, a 
closed form equation that describes the decay of the charge is developed. The 
new equation reduces to the single exponential decay for a material classified as 
a good conductor. For a good insulator the new equation reduces to the known 
"hyperbolic decay law" during the initial decay time. 

Traditional Charge Decay Development 

Most materials can be classified as "simple" materials and this paper deals only 
with such materials which arc defined as isotropic, linear and homogeneous [1]. 
Isotropic means that the polarization vector points in the direction of the elec- 
tric field vector E. Linear means the electric susceptibility is a scalar constant so 
the polarization vector and electric field vector are linearly related, and hence, 
the material's permittivity e ^ e(E) is not a function of the electric field. Ho- 
mogeneous means no spatial gradients and hence Ve = and Va = 0. Within 
a material Gauss' law gives the relationship between the electric field E at a 
point and the net charge density p at that point as [1, 2] 

V • eE = p (1) 

where V is the usual "del" operator ([!]; pp. 39-43). The general continuity 
equation includes source terms (e.g., rate of charge generation) and sinks terms 
(e.g., rate of charge recombination) [3, 4]. At times these rate terms must be 
included in the equation of continuity. However, if sources are equal to sinks 
at all points or if these sources and sinks do not exist, the net charge fiux or 
current density vector must satisfy the simplified (source and sink free) equation 
of continuity [1, 2] 

|^ + V.J = 0. (2) 
A simple material follows Ohm's law given by [1, 2] 

J = ctE. (3) 
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For a simple material that follows Ohm's law, combining (2) and (3) gives 



^ + aV • E + E • Vcr = 0. (4) 
at 

As noted above for a simple material Vcr = and Ve = so (1) and (4) yield 

The traditional method to solve (5) is simply to state when charge density p is 
added that a and e are still constants and define the electrical relaxation time 
of the conductor as 

r. = l (6) 

and integrate (5) from time t = Q when the charge density is po to general time 
t obtaining 

=poexp(-VTc). (7) 

If a unipolar charge is placed in a material, the coulomb force will push the 

charges away from each other. Eventually all charges will arrive at the surface 
so (7) indicates the net charge p at all points inside the material will decay to 
zero in approximately 5rc. The electrical conductivity of copper is 5.8 x 10^ 
S/m and the dielectric constant is 1 giving Tc = 1.53 x 10"^'^ second, whereas 
for fused quartz u = 10~^^ S/m and = 5 giving Tc = 51.2 days; so, in many 
applications copper is considered a good conductor moving any injected charge 
to the surface of the conductor within lO"'^^ second whereas quartz is considered 
an insulator suggesting that for the first several hours any charge introduced in 
quartz will behave as if it had remained essentially wherever it was placed [2] . 

In practice, it is found that metals and other good conductors follow (7) but 
that for insulator materials the initial decay appears hyperbolic with a decay 
rate dependent on the initial charge present. As a result (7) is clearly not valid 
for both conductors and insulators. Later it will be argued that although (5) 
is the proper partial differential equation to describe the relaxation of the net 
charge with time, a more thorough investigation of Ohm's law shows that when 
charge is added to a material a = a{t) and, hence, the traditional method of 
solving (5) described above to obtain (7) is not actually valid. 



Ohm's Law Revisited 

Today there is sufficient evidence to indicate matter is made up of neutral atoms 
and groups of atoms called molecules. Atoms are believed to be made up of a 
nucleus containing protons and neutrons with the nucleus surrounded by elec- 
trons. Neutral atoms have equal numbers of protons and electrons. An electron 
has a negative charge measured at -1.6 x 10~^^ coulomb (C) and a proton has a 
positive charge of +1.6 x 10~^^ C. It is very difficult to remove a proton from an 
atom, but an electron can be more easily removed. As a result, positive charges 
occur in matter when electrons are removed from neutral atoms; and negative 
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charges occur when electrons are added (attached) to neutral atoms or remain 
free. This information and especially the fact that the intrinsic composition of 
a material is made up of neutral atoms and neutral molecules is important in 
understanding the natiue of Ohm's law. 

The general charge flux equation for a single species of free charge is well 
known and given by ([3]; see their Eq. 18) 



Equation (8), which is also referred to as a constitutive law [5] can be obtained 
the traditional way by studying individually each transport process that can 
occur in a volume element at a point and then adding all these processes together 
[6] or by applying the contiguous-collision-averaging method to a point in space 
[7] . The nomenclature used here is that given in [7] in which the sign of charge 
Si is explicit so that the charge density is [7] 



where is the number density of the charge species having a charge of magni- 
tude Qi where Qi = ZiqQ with qq = +1.6 x 10^^^ coulomb being the fundamental 
unit of charge and Zi being the number of fundamental units. The term Si is 
the sign of the charge where Si = 1 for positive charge and Si — -1 for negative 
charge. As a result qiUi is always a positive quantity and the sign of pi is de- 
termined by the sign of Si. The introduction of the term Si allows the mobility 
constant bi to always be introduced as a positive qTiantity where the velocity of 
the charge is given by = SibiE. The conductivity of the i*'* species is [7] 



The diffusion constant of the i charged species is and the thermophoresis 
coefficient is G,, but these need not be considered further here as this paper is 
restricted to homogeneous materials (Vp; = 0) at uniform temperature (VT = 
0). The drift velocity v^o in (8) of the material can often be important in liquids 
and gases and is zero for a stationary solid, liquid or gas. 

Usiially more than one species of charge exits in a material. For example, in 
a neutral gas ionization removes electrons from neutral molecules producing free 
electrons as one species and positive ions as another species. The electrons can 
then attach to neutral molecules and produce negative ions as a third species. A 
similar result occurs in liquids due to dissociation. In solids there can be many 
(conductor), some (semiconductor) or few (insulator) free electrons. Electrons 
freed from atoms in solids also produce positive charged atoms, but these atoms 
cannot physically move. However, by considering these charged atoms as holes, 
when an electron moves from its neutral atom to a previously charged atom 
(hole) the original hole disappears but a new hole is created. This effective hole 
motion is treated as a moving charged species [8]. In general the total charge 
flux is thus given by the sum of all the charge fluxes [7, 9] 



(8) 



Pi = SiQiUi 



(9) 




(10) 




(11) 
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Consequently, (8) and (11) combine to define the total charge flux as 

J = ^ PiVdo + ^ a,E - ^ A Vpi - J2 GiPiVT. (12) 

i i i i 

In (12) the total or net volume charge density in the material at any point is 
given by [7] 

P = Y.Pi- (13) 

i 

Similarly the total electrical conductivity of the material at any point is [7] 

a = (14) 

i 

With (13) and (14) it is possible to write (12) as [7] 

J = pvdo + aE-J2 D^Vp, - (^ G,pi)VT. (15) 

i i 

As pointed out in [7] many materials have no convoctivc drift v^o = 0, are 
homogeneous Vpi = (i.e., are a simple material) and have a uniform tem- 
perature VT = 0. For these conditions (15) reduces to Ohm's law as given by 
(3). Note that nothing in what has so far been presented requires the total 
electrical conductivity (14) to be a constant with time. However, as Maxwell 
pointed out. Ohm's law would have little scientific value if a specific property 
of the material (i.e., the intrinsic conchictivity) could not be defined [10]. This 
intrinsic electrical conductivity, which is not a function of time, will be clearly 
defined in the next section. 

Consequences of Ohm's Law 

In order to understand the concept of introducing free charge into a material 
the concept of intrinsic conductivity must first be understood. The intrinsic 
conductivity is the electrical conductivity of a material before the introduction 
of external free charges. This can be best understood by examining the impli- 
cations of (13) and (14) regarding a material's p and a before external charge 
is introduced (Case 1). Then the complication of introducing a single species of 
free charge (Case 2) can be easily followed. 

Case 1: (Charge neutrality; intrinsic conductivity) 

Because all atoms have equal amounts of positive and negative charge and be- 
cause all materials are made up of atoms, all materials in their equilibrium or 
intrinsic state will be neutral. Therefore, the intrinsic net charge must consist 
of equal numbers of positive and negative charge; consequently, the intrinsic 
net charge must be zero. This analysis can be easily extended to any number 
of sets of different charge species but consider here the simplest situation were 
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only one set of each charge (i.e., two charged species, one positive and the other 
negative) occurs. When no other charges are present except those from the 
material's atoms and molecules, the equilibrium density is from (9) and (13) 

p{eq) ^ p,n^ P+ + P- = s+q+n+ + S-Q-TI- = s+(g+n+ - g_n_) (16) 

where it has been recalled that s+ = — s_ = 1. From (10) and (14) 

(j(eg) am = (t+ + a- = s^g+n+5+ + s^q^n^b- (17) 

However, atomic charge neutrality requires the material to have no net charge 
{pm = ) in equilibrium [8], SO pin — in (16) giving (for the special case 
n+ = n- i.e., q+ = q-) 

q+n+ = q-n- = qmrim (18) 

On the other hand, these free charges give rise to the conductivity of the mate- 
rial, so 

f^m = s\qmnm{b+ + 6_) (19) 

must be the intrinsic electrical conductivity of the material. This same proce- 
dure applies when n+ ^ ri- (i.e., some liquids, gases), but then q+ ^ q-, so 
(19) remains in the form of (17). However, due to the atomic charge neutrality 
principle, in equilibrium = will always be true. 

Case 2: (Perturb with a single charge species pp) 

Next consider what happens when extrinsic charge Pp is uniformly placed in the 
material. For this case (9) and (13) require 

P = Pm+ Pp = Pm + SpqpHp. (20) 

If only interactions (collisions and transport) without reactions (changes in in- 
trinsic ion concentrations) occur between the material and the inserted charge 
then remains fixed. Hence, with the assumption of interaction without re- 
action Pm = 0, so (20) becomes 

p = Pp = SpqpUp. (21) 

Furthermore, (10) and (14) require 

o" = (Jm -I- (Jp = (Jm -h SpqpUpbp = am + SpPpbp. (22) 

Clearly, when the perturbation charge pp is placed in the material, it is expected 
that coulomb repulsion will decay pp away in time. As a result, it is expected 
that Pp = Pp{t) which implies tXp = (Tp{t), and, therefore, a = a{t) in (22). 
The conclusion is clear: when charges are added to a material the electrical 
conductivity of the material is changed; and not until the charges fully decay 
does the electrical conductivity of the material arive back at its intrinsic value. 
As a result, a in (5) is not independent of time; hence, (7) is not the general 
solution to (5). 
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General Charge Decay Equation 

As discussed above, it is clear that (21) and (22) are required in (5); so, (5) 
becomes 

The material time constant can be defined with the intrinsic material conduc- 
tivity as 

Tm = —- (24) 

The constants in the third term of (23) can be lumped together as 

/? = (25) 

With the constants (24) and (25) the partial differential equation (5) i.e., (23) 
becomes 

^ + £^+PpI = 0. (26) 



The solution to (26) is given by [11] 



f^TmPp + 1 (iTmPpn + 



-cxp(-i/T„). (27) 



It is advantageous to define a time constant specified by the initial perturbation 
charge as 

With some rearranging (27) can now be written with the aid of (28) as 

= P.oexp(-Vr^) 

l + ^[l-exp(-VT„)] 

Tp 

which is the correct solution to (5) or (23). Because (29) was derived without 
specifying the state (solid, liquid or gas) and the nature (good conductor, good 
insulator, etc.) of the material, (29) can be applied to all materials from very 
good conductors to very poor conductors (i.e., good insulators). 

Definition of a Good Conductor 

One way to define a good conductor is to require the third term in (26) to be 
negligible compared to the second term. By this definition the criteria for a 
good conductor is Pp/Tm '> 0Pp or 1 ^ ^TmPp- However the largest value is at 
PpQ, so, the criteria for a good conductor must be 1 » (irmPpO- With the aid 
of (28) this criteria becomes Tp ^ Tm- With this good conductor criterion (i.e., 
restriction) (29) reduces to 

Pp(*) ~ Ppoexp(-i/rm), only valid if Tp > r^. (30) 



7 



Equation (30) is similar to (7) except the material's intrinsic conductivity am 
must be specified since, as shown in this paper, a is not a constant. 

Definition of a Good Insulator 

One way to define a good insulator is to require the second term in (26) to be 
negligible compared to the third term. Hence, for a good insulator pp/Tm <C Pp^p 
or 1 <C f^TmPp is required. However, this can only be satisfied at the early stage 
of charge decay, because later, as t oo so must pp decay to zero, and, therefore, 
the requirement 1 <C (ir-mPp eventually always fails. At early time, the only time 
when the requirement can be met, the largest value is at ppo; so, the criteria for 
a good insulator is 1 <C jSTmPpO- With the aid of (28) this criterion for a good 
insulator becomes Tp <C Tm. 

Comparison to Other "Insulator Laws" 

The neumerator and denominator of the single species charge decay equation 
(29) can be multiplied by ex.p{t/Tm) to yield 

Pp{t) = r — — (31) 

(1 H )exp(i/rm) 

Equation (31) has a form somewhat similar to the so-called modified hyperbolic 
law [12] as developed by Vellenga and Klinkenberg [13] in that both equations 
contain an exponential term and a time-constants-ratio term in their denomi- 
nators. However, the two equations are quite different in decay character in the 
range lOrp. However, (29), and hence also (31), reduces to the modified 

hyperbolic law for t„i ^ Tp. 

It is interesting to note that, in general, for t ■^Tm the exponential is given 
by exp(f/Tm) PS 1 + t/Tm which when inserted into (31) gives 

Pp{t) « /''° ; only valid if t < r„. (32) 

For a very good insulator Tp ^ Tm and (32) reduces further to 

pp{t) w — ; only valid if Tp <^Tm, and only if t <^ Tm- (33) 

Tp 

Although (33) is the hyperbolic law of Bustin, et al. [14], the hyperbolic law is 
often (apparently improperly) assumed to be applicable over all time. As noted 
above, (33) shows the hyperbolic law is only valid for t <^ Tm- 

Finally, note that = f3ppo = [{spbp)/e]ppo; and, since ppo = SpQpUp, the 
product in Tp always contains s^. Thus, Tp is always a positive quantity. 
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Discussion 

Equation (29) derived in this paper is a general equation for charge decay of 
unipolar ions which have an initial charge density ppo and which have been uni- 
formly introduced into any material: conductor or insulator; solid, liquid or gas. 
It is assumed that the material media into which the charge is introduced can be 
classified as a simple material; i.e., the material's electrical properties must be 
isotropic, linear and homogeneous. In the derivation of (29) it is assumed that 
the inserted charge pp only interacts with the material through collisions and 
transport and not through reactions which can change the intrinsic ion concen- 
trations of the material. In (29) the intrinsic material electrical relaxation time 
constant t„i is given by (24) where e is the material's permittivity and where 
a^n is the intrinsic electrical conductivity of the material. The time constant 
ratio in (29) is 

7m _ bpSpPpO ^g^^ 

where bp (always a positive quantity) is the electrical mobility of the unipolar 
ions of initial volume charge density Ppo and Sp is the sign of the charge (either 
-1-1 or -1) of the unipolar ions. The initial charge density of the unipolar ions in 
the material is 

PpO = SpQpHpO (35) 

and depends on the number of ions per unit volume Upo, and the magnitude of 
the charge qp of the ions where cjp = ZpQpQ and where qa = 1.6 x 10^^'^ C is the 
fundamental unit of charge and Zp is the number of fundamental units. Note 
that SpPpo is always positive. 

The definition and distinction between a good conductor and a good insu- 
lator can be defined using the the two electrical relaxation times given in (24) 
and (28). Namely, the criterion for a good conductor is Tp ^ and for a 
good insulator is Tp ^ . Since the mobilities are often of the same order of 
magnitude in a given material, these definitions also define a good conductor as 
rim ^ np and a good insulator as Um ^ np. These definitions and (29) are use- 
ful in specifying the type of charge decay (effectively exponential or early time 
hyperbolic). Of course, the basic application definition of a good conductor still 
remains Tm <C Tuse and for a good insulator Tm ^ Tuse where Tuse is the use 
time for the process, investigation, operation, etc. 

Summary 

A general equation (29) has been developed which describes the time decay 
for a uniform distribution of unipolar charge inserted into any simple material 
(conductor or insulator; solid, liquid or gas) with an initial charge density PpQ. 
The relaxation is controlled by both the intrinsic properties (conductivity a^n 
and permittivity e) of the material, and the extrinsic properties of the inserted 
charge, i.e., the initial volume charge density ppo and the mobility bp of the 
inserted charge. In (29) the charge density pp is assumed not to chemically react 
with the material {am invariant to pp) and decays in time t depending on the 
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initial charge density ppo defined by (35), on the material's intrinsic relaxation 
time constant Tm defined by (24), and on the perturbation time constant Tp 
defined by (28). 
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